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i— 1 \ Abstract 

The conformational vibrations of Z-DNA with counterions are studied in framework 

r^M | of phenomenological model developed. The structure of left-handed double helix with 

counterions neutralizing the negatively charged phosphate groups of DNA is considered as 

the ion-phosphate lattice. The frequencies and Raman intensities for the modes of Z-DNA 

with Na + and Mg 2+ ions are calculated, and the low- frequency Raman spectra are built. 

At the spectra range about the frequency 150 cm -1 new mode of ion-phosphate vibrations 

is found, which characterizes the vibrations of Mg 2+ counterions. The results of our 

£> , calculations show that the intensities of Z-DNA modes are sensitive to the concentration 

of magnesium counterions. The obtained results describe well the experimental Raman 

^^ ■ spectra of Z-DNA. 
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: 1 Introduction 

co ; 

DNA is macromolecule with the double helix structure which may adopt different forms as 
a response to changes of environmental conditions [THE]. Under the natural conditions the 
double helix is right-handed due to the stabilization by metal ions (counterions) that neutralize 
the negatively charged phosphate groups of the macromolecule backbone. Increasing the con- 
centration of counterions the double helix may take the left-handed conformation with zigzag 
like backbone (Z-form), which is significant in many biological processes [6j. The study of 
counterion influence on the structure and dynamics of the left-handed form of the double helix 
is important for the understanding the mechanisms of DNA biological functioning. 

The counterions of different type may be localized in different manner with respect to 
the double helix phosphate groups. The experimental data for solid samples of DNA show 
that monovalent counterions are usually localized near the oxygen atoms of the phosphate 
groups from the outside of the macromolecule [HE], while the counterions of higher charge 
may neutralize the phosphate groups of different DNA strands or bind to the nucleic bases 
[9]. In the solution the counterions do not have defined positions, and they form a shell near 
DNA surface [T0H13] . In experiment this shell is observed as a cloud around the double helix 
[HHT7] . The molecular dynamics simulations of DNA counterions in water solution show that 
the residence time of counterions near the phosphate group of DNA backbone is about 1 ns 
[H1[T9]. m this time scale the counterions may form a regular structure around the double 
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helix. Such structure of counterions and DNA phosphate groups may be considered as a lattice 
of ionic type (ion-phosphate lattice). 

The DNA ion-phosphate lattice is expected to have properties of ionic crystals, therefore 
it should be characterized by counterion vibrations with respect to the phosphate groups (ion- 
phosphate vibrations). The ion-phosphate vibrations must be coupled with the internal vibra- 
tions of the double helix, since the dynamics of the ion-phosphate lattice is the part of DNA 
conformational dynamics. The determination of DNA ion-phosphate modes is of paramount 
importance for the understanding of counterion influence on the structure and dynamics of 
DNA double helix. 

The DNA ion-phosphate vibrations may be observed in the low-frequency spectra (IO-t-200 
cm -1 ) where the ion vibrations in ionic crystals and electrolyte solution are prominent [201I2T] . 
In this spectra range the modes of DNA conformational vibrations are also observed [221425] . 
which are characterized by the displacements of the atomic groups in the nucleotide pairs 
(phosphates, nucleosides and nucleic bases) from their equilibrium positions [26H29] . The 
conformational vibrations of the right-handed forms of the double helix are described well within 
the framework of the phenomenological approach [2"oT - |2"9"] . In our previous works [3"0"H3"5] this 
approach has been extended for the case of vibrations of the monovalent couterions with respect 
to the phosphate groups, and the ion-phosphate vibrations are determined for the right-handed 
double helix with Na + , K + , Rb + , and Cs + counterions. The results have been showed that 
the frequencies of ion-phosphate vibrations decrease as counterion mass increases from 180 to 
90 cm -1 . The character of DNA conformational vibrations is found to be very sensitive to the 
counterion type and localization with respect to the double helix atomic groups [Ml 135] . 

In case of Z-DNA the character of conformational vibrations may be essentially different 
because the structure of the left-handed double helix has the unique features. The phosphate 
groups of the left-handed double helix form a backbone with zigzag shape and the nearest base 
pairs have different overlapping. Therefore, the nucleotide pairs form the dimers which are the 
elementary units of Z-DNA [T|. The overlapping of the base pairs inside dimer is stronger 
than between dimers. In the same time, the dimeric structure of Z-DNA is also stabilized 
by counterions (usually Mg 2+ ) that may be localized between phosphate groups of different 
strands of the double helix. The formation of dimers by nucleotide pairs and the localization of 
counterions in the DNA minor grove may be essential for the low-frequency modes of Z-DNA. 
The experimental data show that in the low-frequency spectra of Z-DNA new modes near 150 
and 42 cm -1 are observed [SB]. The origin of these modes is not determined yet. To describe 
the low-frequency spectra of Z-DNA the approach [301435] should be extended considering the 
features of the structure of zigzag like form of the double helix and possible cases of counterion 
localization with respect to the phosphate groups. 

The goal of the present work is to find the modes of optic type for Z-DNA with counterions 
and to describe the conformational dynamics of the left-handed double helix. To solve this 
problem in section 2 the approach for the description of dynamics of DNA ion-phosphate 
lattice [30H35] is extended for the case of Z-DNA. In section 3 the frequencies and the Raman 
intensities of the modes of DNA conformational vibrations are calculated and the low-frequency 
spectrum of Z-DNA is built. Using the calculated spectrum the origin of new modes of Z-DNA 
observed in the experimental spectra is determined. 
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Figure 1: Ion-phosphate lattices of Z-DNA. (a) Zigzag-like shape of Z-DNA backbone, (b) 
Model for conformational vibrations of the left-handed double helix with counterions. The 
monomer links of the model consisted of two nucleotide pairs (dimers) are showed by dotted 
frames. 



2 Model for conformational vibrations of ion-phosphate 
lattice of Z-DNA 

To describe the conformational dynamics of DNA ion-phosphate lattice the typical cases of 
counterion localization with respect to the phosphate groups should be considered. Usually the 
monovalent counterions are localized from the outside of the double helix, where one mono- 
valent counterion neutralizes one phosphate group (single-stranded position of counterion) . In 
the same time, due to the zigzag like shape of the backbone of Z-DNA (Fig. la), the distance 
between phosphate groups of different strands is rather short that makes possible the localiza- 
tion of counterions in the minor grove of the left-handed double helix between phosphates of 
different strands (cross-stranded position of counterion). The cross-stranded position is more 
favorable for bivalent counterions (Fig. la). 

To describe the conformational vibrations of DNA double helix the phenomenological ap- 
proach [26H28] is used. In framework of this approach the phosphate groups (PO4+C5') and 
nuclosides are modelled as masses m® and m, respectively. The nucleosides rotate as the phys- 
ical pendulums with respect to the phosphate groups in plane of nucleotide pair. The physical 
pendulums are characterized by reduced length I. The nucleosides of different chains are paired 
by H-bonds (Fig. lb). The motions of structural elements of the monomer link are considered 
in the plane orthogonal to the helical axis (transverse vibrations). The longitudinal vibrations 
of the macromolecule atomic groups have much higher frequencies and are beyond the scope of 
this work. 

According to the dimeric structure of the zigzag lake double helix the monomer link of 
the model of Z-DNA consists of two nucleotide pairs. The nucleosides in dimer usually have 
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Figure 2: Monomer link of the left-handed double helixes (dimer of nucleotide pairs). I is 
reduced length of pendulum-nucleoside; 6*0 is equilibrium angle; m, m , and m a are masses of 
nucleosides, phosphate groups, and counterions, respectively; X, Y, 9, p, and £ are vibrational 
coordinates of the model (see text). The arrows indicate positive directions of displacements. 



alternated sequence of guanine (G) and cytosine (C) nucleic bases |lj. The displacements of 
nucleosides and phosphate groups in DNA monomer link are described by coordinate Y. The 
coordinate 9 describes the deviations of pendulum-nucleosides from their equilibrium position 
(angle 9q) in the plane of complementary DNA pair. The vibrations of deoxyribose and base 
with respect to each other, inside the nucleoside (intranucleoside mobility), are described by 
changes of pendulum lengths p. The vibrations of counterion in single-stranded positions are 
described by coordinate £. For description of vibrations of a counterion in cross-stranded 
positions the coordinate Y a is used. The vibrational coordinates of the model and the positive 
directions of displacements are showed on Figure 2. 

Within the framework of introduced model of the dynamics of DNA ion-phosphate lattice 
the energy of vibrations of structural elements of the double helix may be written as follows: 



E = J2 {K n + U n + U n . 



n—l , 



where K n and U n are the kinetic and potential energies of the monomer link n of DNA ion- 
phosphate lattice, U n , n -i is the potential energy of interaction along the chain. The kinetic 
energy of the monomer link may be written as follows: 



K n = K 0n + K, 



(2) 



where Kq h and K an are the energies of vibrations of structure elements in nucleotide pairs of 
n-th dimer and counterions, respectively. The potential energy of displacements of the masses 
in monomer link may be written as follows: 



U n = U 0n + U, 



(3) 



where Uq u and U an are the energy of vibrations of the masses in nucleotide pairs of dimer and 
the energy of vibrations counterions, respectively 

The expressions for K§ n and Uq u for Z-DNA may be written in the following form: 

2 2 

K ^=oJ2 Yy M '-A + m v(plij + l l®ni 3 + 2lt 3 9 mj Y nij + 2b ijPnij Y nij )}, (4) 



1 2 2 1 2 

^On = ~ J^Hni + ^2^ijPnij + PiAij 2 )] + ~ ^Wnlj ~ 6 n2 ^ + ^piPnlj ~ Pn2jf + (5) 



+<7j/0nlj ~ Yn2j) ], 

where If, = %o^-; a^ = sm0 Oi j] b^ = cos 9 0i j] n enumerates dimer s of macromolecule; i = 
1,2 enumerates the nucleotide pairs in these dimers; j = 1,2 enumerates the chain of the 
double helix; a, (Tij, and /3y are the force constants describing H-bond stretching in base pairs, 
intranucleoside mobility, and rotation of nucleosides with respect to the backbone chain in base- 
pair plane, respectively; gg, g p , and g y are the force constants describing the interaction of nucleic 
bases and the interaction between phosphate groups of different nucleotide pairs in dimers. The 
variable 5 n i describes stretching of H-bonds in the base pairs (Fig. 2). In the present work it is 
determined analogically to [23121]: 5 ni w l s n6 nil + l si2 0ni2 + Y nil + Y ni2 + bap n ii + b i2 p n i2- 

The kinetic and potential energies of counterion vibrations in the ion-phosphate lattice of 
Z-DNA may be written as follows: 

K an = ^[{Y n2l + L 2l f + (Y nl2 + 4i 2 ) 2 ] + ^fYL (6) 

Uan = |(^ 21 + &,) + y [{Y an ~ Y n22 f + (Y an + F nll ) 2 ] , (7) 

where 7 and 7 a are the force constants for the counterions in single-stranded and cross-stranded 
positions, respectively; m a and M a are the masses of counterions in single-stranded and cross- 
stranded positions, respectively. In equation (6) and (7) the first terms describe the kinetic and 
potential energy of vibrations of counterions in single-stranded positions, while the last terms 
describe the the kinetic and potential energy of vibrations of counterions in cross-stranded 
positions. 

To describe the low-frequency Raman spectra of Z-DNA we intent to find the optic type 
phonons, which are observed in the experimental spectra from 10 to 200 cm^ 1 . Therefore, 
we will consider the modes characterizing the vibrations of DNA atomic groups inside the Z- 
DNA monomers (dimers of nucleotide pairs) which are prominent in the considered spectra 
range. The vibrations of the pair dimers with respect to each other by our estimations occur 
with essentially lower frequencies (< 5 cm -1 ). Since only the long-range optic vibrations are 
observed in the low-frequency Raman spectra we will consider the limited long-wave (k — >-0) 
vibrational modes of the ion-phosphate lattice in Z-DNA. From the point of the theory of the 
lattice vibrations such approximation is the same as the neglecting interaction along the chain 
[37J. So, in the following consideration we neglect the interaction term: U n ^ n -\ ~ 0. Within 
the framework of this approximation the equations of motions may be written as follows: 



i «„ d a^, _ a^ _ au^ 



dt dq n dt dq n dq n dq n 

where q n denotes some vibrational coordinate in the monomer link of model (Fig. 2). The 
equations of motion (jSJ) in explicit form for the model coordinates of Z-DNA double helix are 
showed in Appendix. 

To estimate the frequencies of Z-DNA conformational vibrations the force constants a = 85 
kcal/molA 2 , (3 = 40(46) kcal/mol, a = 43(22) kcal/molA 2 are used the same as in B(A)-DNA 
j26j|27]. The values of /3 and a in case of the left-handed Z-DNA are different for G and C 
nucleosides, and they are taken the same as in A- and f?-forms, respectively. The constants ge 
and g p are equal to 20 kcal/mol, while the constant g y = 8 kcal/molA 2 , [2SJ[27]. The structure 
parameters of the model (I and 9q) for Z-DNA are determined using the X-ray data (pdb code: 
ldcg) [38]. 

The constants of ion-phosphate vibrations have been determined in our previous works 
[SUES], by using the theory of ionic crystals. It has been showed that the constant of ion- 
phosphate vibrations depends mostly on dielectric constant of DNA ion-hydrate shell, the 
Madelung constant of the system describing the electrostatic interaction of the ion with all 
charges of the system, and parameters of the ion (Pauling radius of cation and oxygen atom). 
The dielectric constant has been determined considering the properties of the hydration shell 
of DNA with the counterions of different type. The Madelung constant has been calculated 
from the structure of the double helix. As the result the values of 7 and 7 a constants have been 
calculated for the case of ion-phosphate lattice of DNA with alkali metal counterions in single- 
stranded position and for magnesium counterions in cross-stranded position. The magnesium 
counterions are considered with the hydration shell because the size of hydrated Mg 2+ ion 
corresponds to the distances between phosphate groups of the left-handed double helix. The 
hydration shell of magnesium ion consists of 4 water molecules that are strongly bond with the 
ion J3U]. Thus, the constants of ion-phosphate vibrations have the following values: 7 = 52 
kcal/mol A 2 in case of Na + counterions [31], while in case of Mg 2+ 7 a = 62 kcal/mol A 2 [35J. 
Using such parameters the frequencies of Z-DNA conformational vibrations are estimated by 
the formulae (JTJ - ©. 

3 The low-frequency Raman spectra of Z-DNA 

According to the character of vibrations of structural elements in DNA nucleotide pairs the 
obtained modes may be classified as the modes of ion-phosphate vibrations (Ion), H-bond 
stretching modes (H), modes of intranucleoside vibrations (S), and modes of backbone vi- 
brations (B). To characterize the influence of counterion mobility the frequencies of Z-DNA 
conformational vibrations for the dimers without Mg 2+ counterions are also calculated. The 
frequency values for Z-DNA conformational vibrations are showed in the Table 1. 

The calculations show that the frequencies of ion-phosphate vibrations depend on counterion 
position. In case of Na + counterions in single-stranded position there are two degenerated 
modes (co Ion i and co Ion2 ) with the frequencies at the top of DNA low-frequency spectra (near 
180 cm -1 ). These modes are characteristic for the both right- and left-handed double helix 
forms. The vibrations of hydrated Mg 2+ ions in cross-stranded position {uJianz) have essentially 
lower frequency that is due to the big size and big mass of hydrated Mg 2+ ion (M a = 96 a.u.m.). 
The vibrations of Mg 2+ counterion in cross-stranded position do not influences the vibrations 



Table 1: Frequencies of Z-DNA conformational vibrations (cm -1 ). The results for two cases 
are showed: Z-DNA with Na + and Mg 2+ counterions, and Z-DNA with Na + counterions only. 
The calculated frequencies are compared with the experimental data for Z-DNA [23~||3"6] and 
with the calculation data for i?-DNA [311. 
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of Na + counterions outside of the double helix. 

The H-bond stretching in nucleotide pairs of Z-DNA are characterized by 4 modes (oohi, 
&H2, &H3, and cuf/4), while in the case of the right-handed form of the double helix there are 
only two modes. Additional two modes of H-bond stretching ujh2 and ujh3 appear due to the 
formation of dimers of nucleotide pairs in Z-DNA structure. The modes ujhi and uh2 have 
practically the same frequency for B- and Z-DNA. The modes ojh3 and u)H4 of Z-DNA are 
sensitive to Mg 2+ counterion in cross-stranded position. 

In case of left-handed double helix there are 3 modes of intranucleoside vibrations ojgi, ^S2, 
and Us3 with rather close frequency values (about 45 cm -1 ). The modes of intranucleoside 
vibrations depend essentially on the double helix form. In case of 5-DNA there is only one 
mode of intranucleoside vibrations usi, which has essentially higher frequency than respective 
modes of Z-DNA. The mode Us3 shifts to lower frequencies if there is no Mg 2+ counterion in 
cross-stranded position. 

The modes of backbone vibrations (u>bi, ^B2, ^B3, and coba) are the lowest in DNA low- 
frequency spectra. In case of the ion-phosphate lattice of Z-DNA there are two modes (ujbi, 
ub 2 ) near 25 cm -1 and two modes (cjb3, ^ba) with close frequency values (about 10 cm -1 ). In 
case of the ion-phosphate lattice of i?-DNA there are only two modes ujbi and u>B2 with close 
frequency values. This difference of B- and Z-DNA spectra is the manifestation of the dimeric 
structure of the left-handed structure of DNA double helix. 

To compare our calculations with the experimental spectra the Raman intensities are calcu- 
lated using the approach developed in [3U|35]. As the result the low- frequency Raman spectra 
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Figure 3: The low-frequency Raman spectra of Z-DNA with different filling of the ion-phosphate 
lattice by Mg 2+ counterions. a) Backbone vibration range (0-=-40) cm _1 (B). b) Intranucleoside 
vibrations range (40-=-60) cm -1 (S). c) H-bond stretching vibrations range (60-=-130) cm" 1 (H). 
d) Ion-phosphate vibrations range (130-=-200) cm -1 (Ion). 



of DNA in Z-form of the double helix are built for different degree of filling (P) of the ion- 
phosphate lattice by Mg 2+ counterions (Fig. 3). The parameter P equals to 1 if all dimers of 
Z-DNA contain Mg 2+ counterions, and P equals to if there is no counterions in dimers. The 
mode intensities are normalized per intensity of the mode ujhi- The halfwidth of spectra lines 
is considered equaled to 5 cm -1 

The obtained spectrum is divided into four ranges with similar intensities of the modes 
of DNA conformational vibrations (Fig. 3a - 3d). The range (0-=-40) cm -1 (Fig. 3a) is 
characterized by the modes of backbone vibrations (B). This modes are the most intensive 
in DNA low-frequency spectra because their energy is lower than k^T at room temperatures. 
Due to the high intensity the modes of backbone vibrations are observed in the experimental 
spectra as a peak near 25 cm -1 [32] • Increasing the filling of the ion-phosphate lattice by Mg 2+ 
counterions in the dimers the intensity of this band decreases. 

At frequency range (40-r60) cm -1 (Fig. 3b) there is a band characterizing the intranucleoside 
vibrations (S). The intensity of this band is essentially lower than the intensity of B band 
formed by modes of backbone vibrations. Therefore in the experimental spectra [36J S band 
is observed as a weak mode near 45 cm -1 . The dependence of intensity of this band on filling 
parameter P is weak. 

The frequency range (60-=-130) cm" 1 (Fig. 3c) is characterized by two bands 70 and 110 
cm -1 that are caused by the modes of H-bond stretching in the base pairs (H). The intensity 
of the band 70 cm -1 decreases as the filling of the ion-phosphate lattice by Mg 2+ counterions 
increases, while the intensity of the band 110 cm -1 remain the same. The comparison of the 
calculated spectra (Fig. 3c) with the experimental Raman spectra [23,35] show that at these 
frequencies there are two strong modes. The frequency and intensity of the observed modes 
agree with our calculations. 

At the spectra range (130-=-200) cm -1 (Fig. 3d) there are two bands 150 and 180 cm -1 
characterizing the ion-phosphate vibrations of Mg 2+ and Na + counterions, respectively (Ion). 
The dependence of intensity on filling parameter P is essential only in case of the mode of 
ion-phosphate vibrations of Mg 2+ counterions. The comparison with the experimental spectra 
show that in the Raman spectra of Z-DNA very weak mode near 153 cm -1 is observed [36J 
that corresponds to the calculated band 150 cm" 1 . The vibrations of Na + counterions are not 
observed in the low-frequency spectra of Z-DNA. However, in the spectra of dry DNA films 
the mode resembling the mode of ion-phosphate vibrations has been detected [2"2" | |2"5 ] I3"0~] , I3T]. 

4 Conclusions 

The conformational vibrations of Z-DNA are studied considering the double helixes with coun- 
terions as the lattice of ionic type (ion-phosphate lattice). To find vibrational modes of the 
ion-phosphate lattice the phenomenological model is developed taking into consideration the 
feature of Z-DNA structure. In the model the monovalent counterions are localized outside of 
the double helix, and bivalent counterions are localized between phosphate groups of nucleotide 
dimers of Z-DNA. Using the developed model the frequencies and the Raman intensities of vi- 
brational modes are calculated for DNA with Na + and Mg 2+ counterions. The results show that 
the obtained frequencies of H-bond stretching in nucleotide pairs and the frequencies of intranu- 
cleoside vibrations differ from respective frequencies in case of the right-handed DNA double 
helix. The reason of such difference is in the dimeric structure of the left-handed Z-DNA. New 
mode of ion-phosphate vibrations about the frequency 150 cm -1 is determined, which charac- 



terizes the vibrations of Mg 2+ counterions. The intensity of this mode is rather large, and it 
is observed in the experimental Raman spectra of Z-DNA [36J. The modes of ion-phosphate 
vibrations in case of Na + counterions (180 cm -1 ) are weak and do not observed experimentally. 
The modes of Z-DNA conformational vibrations are very sensitive to the concentration of Mg 2+ 
counterions in the solution. The developed model describes the conformational vibrations of 
Z-DNA with Na + and Mg 2+ counterions and allows determine the origin of new modes in the 
experimental low-frequency Raman spectra of Z-DNA. 

The present work was partially supported by the Grant for the Young Scientists of the NAS of Ukraine 
0112U005857 and by the project of the NAS of Ukraine 0110U007540 

Appendix: Equations of motion 

For a more convenient form of the equations of motion the following variables are used: Y™ = 

~yn I ~yn „,n \rn x/~n ()n cm _i_ nn „n on _ nn n n n n _i_ n n r n 

z ll "T" l 22i y C — J ll 2 22; u c — ^11 "T" t7 22; '/c — "ll u 22i Pc ~ Pll "•" P22i ' c ~ 

n n n n \z~n \rn _i_ \sn „.n \sn yn nn nn i nn „n nn nn „n 

Pll ^225 x g — 1 21 "r 1 125 iJg ~ 1 21 1 125 °g ~ v 2\ ' t7 12? 'Ig ~ V 2\ v \2i Pg ~ 

P21 + Pi2> r g = P21 ~ P121 fi = £21 + tv2i £2 = £21 ~ £12- As the result the equations of 
motion (JSl) for the optic long wave vibrations are split into two subsystems of coupled equations: 



Yc n + ^P n c + =g^ + jfcl^? ~ 9yc(Y c n - r;) + 

a c (Y; + Y c n + IT + ITc + hP n g + bcP?) = ; 

Vn 1 ™>gb g ^ n 1 m-glg qu , m a l\rn _i_ cn\ _ (Vn Vn\j_ 
X 9 "•" ~Wg~P9 "•" ~Mi a 9 ~T WgVg ^^l) dygK^c r g) + 

a g (17 + Y c n + V g 9 n g + l' c &2 + b g p n g + b cP n c ) = ; 

Q n + yn^ + ^n _ ^nn _ fln) + Q ^ (yn + yn + ^n + ^n + ^n + ^n) = Q . 
§n + yna t + m _ ^Qn _ ty + a< *£ (yn + yn + qgn + ^n + bgf) n + ^) = Q . 

Pi + 6,17 + a 9 p n g - g eg (p n c - ft) + a g ^ (Y? + Y? + l s g 9 n g + l'& + b g p n g + b c p n c ) = ; 
Pl + b c Y c " + a cP : - g pc (p n c - p n g ) + a c ^t {?? + Y? + Wg + Wc + hp n g + b c p?) = ; 

yn + fflij* + ^n _ H^ {2y n _ y?) + ^ ( y » + y „) + 

« c {y n g + y n c + l s g Vg + IcVc + b g r n g + 6 c r») = ; 

« 9 (l/ 9 n + V n c + l s 9 V n 9 + ItVc + b 9 r r ; + b c r:) = ; 

Vg + vTtg + « + ^ (»?? + O + « 9 ^7 (^ + y U o + W + l o< + ^ r ," + ^ r ") = ; 
ff c + yTt + PcVc + 98c {€ + Vg) + a c %± [y n g + y n c + l^ + ^c + K^ + b c r^) = ; 
r n g + b g y n g + ° a r n g + 9 Q g (r n c + r») + a g ^ (y n g + y n c + Iffi + l 8 c tc + h^g + b c r?) = ; 
r n c + b c y n c + o c r n c + g pc (r™ + r«) + a c ^ (y n + y n c + l s g rf g + l s c rf c + b g r n g + b c rf) = ; 

2y: + 2 la0 (2Y a n -y?)=0, 

where a g = a/M g , a c = a/M c , (3 g = (3/m g l 2 g , (3 C = P/m c l 2 , a g = a/m g , a c = cr/m c , 70 = 7/ m a, 
7ao = la/M a , g eg = ge/m g l 2 g , g 6c = ge/mj 2 c , g pg = g p /m g , g pc = g p /m c , g yg = g y /m g , g yc = 
g y /m c , a c = sm9 0c , a g = sm9 0g , b c = cos9 0c , b g = cos9 0g , l s c = l c sm9 0c , l s g = l g sin9 0g . 
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The first subsystem of equations of motion describes symmetric vibrations of nucleotides 
in the dimers, while the second subsystem describes antisymmetric vibrations of nucleotides. 
The obtained equations of motions are solved using the substitution: q n = q n exp(iut), where 
q n and uj are amplitude and frequency for some coordinate of vibrations, respectively. As the 
result the equation for frequencies of long-wave vibrations of Z-DNA are obtained. From the 
first subsystem of equations of motion the equation for frequencies u>hi, Wff3, ^52, ^S3, ^B2, 
u)Bi, and ouioni is derived. The second subsystem of equations of motion gives the equation for 
frequencies ojh2, ^ha, ^si, ^>bi, ^B3, w/ on2 , and u>i on3 . Using the parameters discussed in the 
section 2 the frequency values of the modes of Z-DNA conformational vibrations are calculated 
numerically. The values of the calculated frequencies are showed in the Table 1 . 
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